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We discuss the procedure to resum large logarithms to all orders in the differential 
distributions for DIS event shape variables —^F-- We describe results for two 
variants of the thrust variable, both denned wrt the boson axis in Breit frame but 
with differing normalisations. 

1 Introduction 

Event shape variables in e + e~ reactions have long been the subject of much 
attention from theorists and experimentalists proving popular tools for the ex- 
traction of a s as well as for the study of non-perturbative (power) corrections. 
In the last few years such attention has also been focussed on similar variables 
in DIS with experimental studies being carried out by both the HI El and Zeus 
collaborations. □ In this regard one area of investigation involves the study 
of differential distributions in shape variables r which have the perturbative 
expansion 

--J- = Ai(r,x)a s + A 2 {T,x)a 2 s + --- (1) 
a ar 

While, in the above equation, the coefficient Ai is obtainable through an 
analytical calculation, A 2 can he provided by NLO Monte Carlo programs 
DISENTS and DISASTER++. @ The variable x is the standard Bjorken vari- 
able. At present, experimental studies are confined to comparing the above 
truncated expansion to data. El The findings are that while such a comparison 
is viable for larger values of r it gets progressively meaningless as we go to the 
small t region. The perturbative results appear to diverge at small t instead 
of turning over as indicated by the data. 

This problem is familiar from e + e~ distributions. Essentially the pertur- 
bative coefficients above have a leading behaviour (in the small r region) 

A n (r,x) ~ Ifr 2 "- 1 (jj +■■■ (2) 

At small r the smallness of a s is more than compensated by the large ln(l/r) 
terms which accounts for the divergent behaviour of fixed order results in that 
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regime. In what follows we shall also refer to the integrated shape cross section 




for which for every power of a s there are up to two powers of ln(l/r) (hence- 
forth denoted by L). The differential distribution can be obtained from R by 
straightforward differentiation. 

The leading (a s L 2 ) n behaviour in the perturbative prediction for R is 
transparently associated with n soft gluon emission and the absence of com- 
plete Bloch-Nordsieck cancellations. Such terms are free from any x depen- 
dence since soft gluon emission does not change significantly the momentum 
fraction of the incoming projectile. However there are other sources of large 
logarithms apart from purely soft emission. These include running coupling 
effects and hard collinear emissions. It turns out that hard collinear emissions 
on the incoming leg in DIS are responsible for terms starting at (a s L) n (single 
logarithmic level)with coefficents which are x dependent. 

Clearly in order to obtain meaningful perturbative results, one is required 
to try and resum these large logarithms to all orders in perturbation theory. 
In this regard the only new feature of the DIS resummation in comparison to 
e + e~ shape variable resummation is the above mentioned x dependent single 
logs. Their resummation is possible since they are of leading log DGLAP type 
and lead to a change of scale in the parton distribution, or the emergence of 
the factor q<, ?'' r ^ 2 \ in the final result, u The other aspects of the resummation 

are mainly familiar from experience with e + e~ variables, though the details 
differ from variable to variable. Note that in writing the factor q(x,rQ 2 ) we 
have also implicitly included the contribution from incoming gluons. 

2 Definitions 

We consider below the following two definitions of the thrust variable in the 
Breit current hemisphere: 



where the unit vector n denotes the photon direction in the Breit frame and the 
sum is over all particles in the current hemisphere H c . The definitions differ 




(3) 



te — 1 - T E 



1 




(4) 



2 



only in normalisation but we find that this significantly affects the resummation 
both in procedure and results. □ The latter definition, involving normalisation 
to the current hemisphere energy, is only infrared safe provided one imposes a 
minimum energy cut-off in H c . 



3 Results and Conclusions 

The resummed result for the contribution to F 2 from events with 1 — T < r <C 1 
(where r is now used to denote tq or te ) is given by 



R(x, Q 2 ,t) = x 



^e 2 q q{x,TQ 2 )+C q {x)+C g {z 



Q.Q 



X(a s ,L) (5) 



Here C q and C g are 0(a s ) constant pieces from incoming quark and gluon 
sectors respectively and they are different for tq and te ■ Additionally one has 
the form factor 



E(a s ,L) = exp[Lgi(a s L) + g 2 (a s L) + a s g 3 (a s L) H ] (6) 

in which we control all terms of the type a™ L n+i Land a"L™. In other words 
the functions 171 and g 2 are explicitly computed □ while 53 is unknown. All 
these functions also depend on the variable under consideration and hence are 
different for the different thrust definitions presentedu As discussed previously, 
all x dependent logarithms have been resummed in the parton density function 
q(x,TQ 2 ) and the resulting form factor is now x independent. 

One can expand this formula to 0(a 2 s ) and compare the results at small 
t to those from fixed order Monte Carlo programs. On doing so we find 
reasonable agreement with DISASTER++ but disagreement _with DISENT. 
The exact nature of this disagreement is detailed in our paperLl 

Our result valid for small r can then be matched to NLO estimates from 
the Monte Carlo programs. The matching procedure extends the range of 
applicability of the resummed calculation to larger r values where non loga- 
rithmic pieces in the fixed order result are significant. The essential idea of 
matching is simple: one simply adds the resummed and NLO calculations and 
then subtracts the pieces corresponding to double counting. Schematically we 
can write 

p _ p , /pNLO pNLO\ ry\ 
-ttmat ^ -flics + [ft — K Tes ) (I) 

In the above i? ma t denotes the matched shape cross section and i? rcs the re- 
summed calculation, while the piece (i? NLO — R^°) takes care of adding terms 
in the NLO result for R that are absent in the resummed result expanded to 
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the same order. In practice matching is a technically involved procedure and 
different matching schemes can be used which differ on how to treat subleading 
logarithmic terms. Our prefered matching prescription is presented in the full 
papertl 

Once the resummation is carried out and matched to fixed order, we have 
the best available perturbative prediction at hand. This prediction is at least as 
good as the NLO Monte Carlo estimate over the entire range of shape variable 
values and is far superior at small r. However we have not addressed another 
potential source of uncertainty, namely power behaved corrections. 

The main effect of power corrections will be to simply shift the perturbative 
distributions by an amount proportional to 1/Q towards larger r values. An 
exception to this rule is the current jet broadening variable where the pertur- 
bative prediction is squeezed and shifted due to non-perturbative effects. The 
amount of the shift for the thrust variables is identical to the ppwer correction 
to the corresponding mean values which have been computed^ 

Other variables that are being studied at the moment include the thrust 
defined wrt the actual thrust axis (normalised to E), the current jet-mass 
and the C parameterD When resummed results for these distributions become 
available they should also be matched to NLO and adjusted for power cor- 
rections after which detailed phcnomcnological analysis and comparisons with 
data should become possible. 
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